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We  study  the  flow  of  a  heavy,  viscous,  possibly  non-Newtonian  axisymmetric  jet  of 
liquid  of  density  p  falling  under  gravity  g  into  a  lighter  liquid  of  density  p.  If  the 
change  in  the  momentum  of  the  entrained  lighter  liquid  is  neglected  the  jet  will 
ultimately  reach  a  modified  Torricelli  limit  with  a  speed  given  by 

U(x)  =  [2^0*]* 

and  an  asymptotic  radius  r  -u 

°(*,=[(dbJ-  1,1 


where  x  is  the  downstream  distance,  dp  =  p—p  >  0  and  2 Q  is  the  volume  flow.  An 
exact  asymptotic  solution  perturbing  the  Torricelli  limit  with  effects  of  surface 
tension,  viscosity  and  elasticity  is  given  in  powers  of  x~i.  An  extended  unsteady 
problem  including  effects  of  entrainment  is  formulated  in  terms  of  nonlinear  ordinary 
differential  equations  which  also  account  for  weak  radial  variations  of  the  velocity 
across  the  cross-section  of  the  jet.  These  equations  are  solved  in  a  boundary-layer 
approximation  which  gives 

|2) 

where  p  is  the  viscosity  of  the  ambient  fluid.  Equation  (1)  is  in  agreement  with 
'  experimental  observations  of  jets  of  liquid  into  air.  Equation  (2)  is  in  agreement  with 
|  experimental  observations  of  jets  of  liquids  into  liquids. 


[l.  Introduction 

‘  A  heavy,  possibly  non-Newtonian,  liquid  of  density  p  is  extruded  from  a  pipe  or 
vessel  of  radius  a0  into  a  lighter  fluid  of  density  p.  It  is  assumed  that  flow  in  the  jet 
axisymmetric  with  velocity 

■  U(z,r)  =  txU(x,r)+trW(x,r) 

and  radiuB  r  *  o(z),  where  x  increases  downward  in  the  direction  of  gravity  (see  figure 
i).  It  is  further  assumed  that  the  viscosity  of  the  jet  is  very  much  larger  than  the 
viscosity  of  the  ambient  fluid.  This  implies  that  the  variation  of  U  is  determined  by 
the  internal  dynamics  of  the  jet,  and  ultimately,  for  large  x,  the  accelerating  jet  will 
thin  out  and  U(x,  r)  will  become  independent  of  r.  The  falling  jet  will  entrain  ambient 
fluid,  however,  even  if  the  viscosity  p  of  the  ambient  liquid  is  small,  and  the  change 
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of  momentum  of  entrained  fluid  will  be  very  important  if  the  density  of  the  entrained 
liquid  is  not  too  small.  Thus  there  are  two  cases: 

(i)  the  entrainment  effect  is  not  important  -  this  case  is  realized  perfectly  for  jets 
into  vacuum,  and  it  is  realized  in  practice  in  the  case  of  jets  of  liquid  into  air; 

(ii)  the  entrainment  effect  is  important  -  this  case  is  realized  when  liquids  are 
extruded  into  liquids. 

The  two  cases  are  qualitatively  different  in  that  the  leading  dynamical  balances, 
in  asymptotic  regions,  are  between  the  relative  weight  of  the  jet  per  unit  length  and 
its  change  of  momentum  in  case  (i),  and  between  the  relative  weight  of  the  jet  and 
the  change  of  momentum  of  the  entrained  liquid  in  case  (ii).  Equations  (1.1)  and  (1.2) 
below  show  that  there  is  a  huge  difference  between  these  two  cases.  Our  experiments 
show  that  both  cases  are  realizable. 

We  treat  case  (i)  under  the  mathematical  assumption  that  the  ambient  fluid  is 
inviscid  and  dynamically  passive.  This  means  that  the  stress  exerted  on  the  heavy 
liquid  by  the  light  one  is  hydrostatic.  The  effect  of  the  hydrostatic  pressure  is  to  reduce 
gravity  g  to  ( Sp/p)g ,  where  Sp  is  the  difference  in  density.  Thus,  any  result  for  flow 
of  a  liquid  in  air  or  a  vacuum  which  depends  on  the  gravity  g  also  describes  the  flow 
in  another  (dynamically  inactive)  liquid  when  g  is  replaced  with  Spg/p.  In  the  case 
of  vertical  jets  into  liquid  we  get  the  effect  of  gravity  as  a  constant  acceleration  of 
the  jet  which  reaches  finally  a  terminal  speed  given  by  the  asserted  modification 
[2(8p/p)gx$  of  Torricelli's  formula  (2^x)l.  The  radius  of  the  jet  must  then  get  smaller 
and  smaller,  reaching  ultimately  the  asymptotic  radius 


We  find  a  formally  exact  asymptotic  solution  of  viscous,  non-Newtonian  equations 
in  powers  of  x~i.  We  work  with  a  nonlinear  jet-shape  equation  derived  from  the 
equation  of  motion.  The  same  results  follow  from  direct  exact  analysis  of  the 
governing  partial  differential  equation.  The  algorithm  for  computing  such  an  exact 
solution  in  the  case  of  Newtonian  jets  into  air  was  indicated  by  Kaye  &  Vale  (1969), 
but  they  obtained  only  the  first  two  terms  of  the  asymptotic  series.  Kaye  &  Vale  found 
that  their  asymptotic  solution  was  good  for  an  oil  of  low  but  not  of  high  viscosity. 
For  the  high-viscosity  oil  they  studied  a  nonlinear  jet-shape  equation  (see  (6.23)) 
which  is  a  special  case  of  (5.18)  and  got  good  agreements.  They  called  the  asymptotic 
method  (II)  and  the  method  of  integration  of  the  jet-shape  equation  (I),  and 
maintained  that  there  was  a  strong  distinction  between  the  two  methods.  In  fact, 
their  method  (I)  is  the  same  as  (II)  when  the  surface-tension  terms  are  added  to  (I) 
and  x  is  large.  Clarke  (1969)  treats  viscous  fluids  with  surface  tension  and  he  has 
derived  the  correct  three-term  asymptotic  expansion  (5.19).  He  uses  a  hodograph 
method  with  streamline  coordinates,  so  his  analysis,  which  is  based  on  some 
unnecessary  assumption  about  small  parameters,  and  the  expression  of  his  results  are 
less  direct  and  appear  to  be  less  rigorous  than  our  analysis  shows  they  actually  are. 

The  fluid  dynamics  of  the  jet  near  the  exit  (say,  within  2  or  3  diameters)  is  more 
complicated  because  the  flow  must  adjust  from  a  distribution  compatible  with  flow 
in  a  pipe  to  flow  in  a  jet.  This  problem  is  usually  studied  with  g  =  0  (see  e.g.  Trogdon 
&  Joseph  1980,  1981).  In  this  case  the  flow  reaches  a  uniform  velocity  and  a  final 
swelled  radius  (larger  than  the  pipe  radius  -  much  larger  for  some  viscoelastic  jets 
(see  figures  2e-g).  So  the  limit  Sp-*  0  might  be  expected  to  coincide  with  the  problem 
of  die  swell  as  it  is  usually  formulated.  There  is  nothing  to  accelerate  the  jet  when 
the  densities  are  matched  and,  in  principle,  the  jet  diameter  will  not  vary  downstream 
(actually  the  jet  will  eventually  be  destroyed  by  instabilities  not  considered  here). 
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It  is  not  uncommon  to  simulate  g  =  0  in  experiments  on  die  swell  by  density 
matching,  but  the  idea  behind  that  experimental  technique  can  be  erroneous  if  it 
neglects  the  effects  of  entrainment  of  the  ambient  fluid  by  viscous  action.  In  fact  the 
rate  of  change  of  the  momentum  of  the  entrained  fluid,  which  in  any  case  is  small 
if  the  density  of  the  entrained  fluid  is  small,  can  be  large  when  density  of  the  lighter 
liquid  is  comparable  to  the  density  of  the  liquid  in  the  jet.  In  this  case  (2)  the  overall 
change  of  momentum  produces  smaller  velocities  and  larger  jets  so  that  the  swelled 
radius  computed  in  case  (1)  for  small  values  of  Spg  may  be  much  smaller  than  those 
that  are  actually  observed.  It  could  be  said  that  there  is  a  ‘die  swell’  which  is 
associated  with  entrainment  of  the  ambient  liquid. 

The  dynamical  effects  of  entrainment  of  ambient  liquid  by  the  jet  are  studied  in 
§§6  and  7.  We  derive  some  nonlinear  ordinary  differential  equations  for  jet  problems 
which  are  axisymmetric  but  unsteady.  These  equations  differ  from  approximate 
equations  used  in  the  problem  of  fibre  spinning,  and  other  thin-fluid  filament 
problems  in  that  the  dynamics  of  the  ambient  liquid  is  represented  and  the  equations 
allow  at  least  for  weak  radial  variations  of  the  velocity  and  pressure  at  each  axial 
station.  In  §7  we  treat  the  problem  of  entrainment  by  approximate  methods  in  which 
the  dynamics  of  the  entrained  liquid  is  presumed  to  be  governed  by  boundary -layer 
equations,  and  the  boundary-layer  equations  are  studied  by  the  method  of  Karman 
and  Pohlhausen.  We  find  that  the  boundary-layer  thickness  increases  like  1/a*  and 
that  ultimately  the  jet  radius  is  nearly  constant  and  is  given  by 


a(x)  =  1171 


\x)  (Spgp  ‘ 


(1.2) 


The  results  from  theoretical  analysis  seem  to  be  in  good  agreement  with  observations 
reported  in  §2  and  discussed  in  §§5  and  8. 


2.  Experiments 

In  our  experiments  we  extruded  heavy  liquids  from  a  round  pipe  of  radius  r  =  o„ 
with  vertical  x-axis  parallel  to  gravity,  into  lighter  fluids  of  different  densities  p.  The 
experiments  were  originally  undertaken  to  study  the  problem  of  ‘die  swell’  in 
non-Newtonian  fluids.  This  led  us  to  certain  problems  in  understanding  the  influence 
of  the  ambient  liquid  on  the  dynamics  of  the  jet.  A  partial  resolution  of  these  problems 
is  established  in  the  theory  developed  in  the  sequel.  To  motivate  our  theory  it  is 
perhaps  best  to  report  some  details  of  the  actual  motions  that  we  observed  in  our 
experiments. 

A  sketch  of  a  typical  experiment  is  shown  below  in  figure  1 .  In  this  work  the  density 
difference  p—p  >  0  is  of  crucial  importance.  We  have  normalized  this  density 
difference  with  the  density  p  of  the  extruded  liquid 


In  a  vacuum  p  ■>  0  and  e  »  1,  e  *  1  when  p  is  for  air  and  e  is  very  small  when  the 
density  p  of  the  lighter  ambient  liquid  is  close  to  p.  Unfortunately  it  is  not  possible 
to  find  ambient  fluids  of  intermediate  density.  We  can  extrude  into  a  gas  or  a  liquid. 
There  is  nothing  in  between.  Fortunately,  we  can  at  least  find  liquids  of  different 
density  so  that  different  but  small  values  e  can  be  used  to  test  our  theory. 

We  worked  with  jets  of  siijoone  oil  and  jets  of  an  elvacite  solution  into  air  Mid  into 
liquids  of  different  density.  The  silicone  oil  is  a  Newtonian  liquid  which  swells  a  little 
at  the  exit  and  the  elvacite  solution  is  a  viscoelastic  liquid  which  swells  a  lot  at  the 


15-S 


446 


D.  D.  Joseph,  K.  Nguyen  and  J.  E.  Malta 


Figure  1 .  Sketch  of  the  falling  jet.  For  large  x  the  jet  radius  and 
the  radial  velocity  tend  to  aero. 


Density 

Viscosity 

difference 

of  ambient 

silicone  1 

fluid 

Ambient  fluid 

(g/cm*) 

(cP) 

Air 

0-943 

— 

100%  alcohol 

0-101 

1-2 

71%  alcohol/ water 

0-068 

2-3 

50%  alcohol/water 

0023 

2-8 

40%  alcohol/water 

0005±0002 

2-9 

Table  1.  Densities  and  viscosities  of  the  ambient  fluid  for  jets  of  silicone' oil  1 


exit.  The  elvacite  solution  was  a  9-8  %  (by  weight)  solution  of  elvacite  in  diethyl 
malonate  (DEM)  and  the  elvacite  itself  is  a  polymethyl  methacrylate  polymer  with 
a  viscosity-average  molecular  weight  of  about  4  x  10*.  First  we  shall  list  the  values 
of  the  relevant  physical  parameters  for  the  liquids  used  in  the  experiments: 
for  silicone  oil  1 

ft  =  10  P,  p  **  0943  g/cm*,  <r  *  21 '2  dyn/em  (into  air); 
for  silicone  oil  2 

ft  «  126  P,  p  «■  0*972  g/cm*,  or  *  21*6  dyn/cm  (into  air). 

The  silicone  oil  was  injected  into  air  and  into  various  mixtures  of  water  and  alcohol. 
The  viscosity  of  these  mixtures  was  of  the  same  order  as  the  viscosity  of  water;  that 
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Density 

Viscosity 

difference 

of  ambient 

silicone  2 

fluid 

Ambient  fluid 

(g/cm*) 

(cP) 

Air 

0-972 

— 

100%  alcohol 

0-187 

1-2 

66%  alcohol/water 

0-102 

2.3 

42%  alcohol/water 

005 

2-8 

35%  alcohol/ water 

0023 

2-8 

27%  alcohol/ water 

0014 

2-5 

Table  2.  Densities  and  viscosities  of  ambient  fluids  for  silicone  oil  2.  The  density 
of  the  alcohol  mixtures  was  determined  by  weighing  of  known  volumes. 


Density  Viscosity 

difference  of  ambient 
elvacite/DEM  fluid 
Ambient  fluid  (g/cm*)  (cP) 

Air  1-067  — 

Water  0-067  1 

25%  glycerol/water  0-005  2 


Table  3.  Densities  and  viscosities  of  ambient  fluids  for  jets  of  elvacite 


is,  many  orders  of  magnitude  smaller  than  the  silicone  oil.  The  density  difference  and 
viscosity  for  the  various  mixtures  used  in  the  experiment  are  given  in  tables  1 
and  2. 

For  the  elvacite  solution 

H  (zero-shear  viscosity)  *  9  P,  p  =  1-067  g/cm®,  cr  =  30-2  dyn/cm  (into  air). 

The  elvacite  was  injected  into  air  and  into  various  mixtures  of  water  and  glycerol. 
The  viscosity  of  these  mixtures  is  also  of  the  order  10~*  P.  The  density  difference  for 
the  various  mixtures  used  in  the  experiment  are  given  in  table  3. 

We  were  unsuccessful  in  finding  values  for  the  interfacial  tension  of  the  silicone 
oil  into  liquid.  Ring-tensiometer  measurements  fail  in  liquid-liquid  cases  in  which 
the  surface  tension  is  too  small  to  measure. 

Some  typical  raw  data  for  the  experiments  are  the  photographs  reproduced  in 
figures  2{a-d)  (silicone  1)  and  2  (e-g)  (elvacite/DEM).  The  mass-flow  rate  for  all  the 
silicone  1  experiments  wasO- 16  g/s(Q  =  l-6/2rrp  *  0027  cm*/s)andforelvacite/DEM 
0-61  g/s  (Q  =  0091  cm*/s).  A  1-78  mm  inner  diameter  nozzle  was  used.  The  axial 
magnification  of  the  true  shape  is  different  from  the  radial  one  owing  to  the  curvature 
of  the  cylindrical  container.  The  jet  shapes  are  also  plotted  in  figure  3  (silicone  graph 
labelled  [2],  [3],  [4],  [5],  [6])  and  figure  4  (elvacite/DEM  graph  labelled  [2],  [3],  [4]). 

The  interpretation  of  the  experiments  is  best  deferred  until  we  have  established 
some  theory.  Direct  comparison  with  theory  can  be  found  at  the  ends  of  §§5  and  8. 
However,  at  this  point  it  is  perhaps  useful  to  remind  the  reader  that  we  are  interested 
in  the  dynamics  of  the  jet  away  from  the  pipe  exit. 

In  the  absence  of  gravity,  e  «  0  the  jet  attains  a  final  uniform  velocity  V  =  Uf  and 
diameter  o  —  a,  (see  e.g.  Reddy  k  Tanner  1978;  Trogdon  k  Joseph  1980,  1981).  In 
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( B ) 

Figure  2.  Photographs  of  falling  jets:  ( a)~(d )  silicone  1 ;  (e)-(p)  elvacite/DEM. 


Figure  3.  Jet  shape  for  silicone  oil  1  injected  into  air  and  various  alcohol  solutions:  [1]  equation 
(4.6);  [2)  ep  =  Sp  =  0-943;  [3)  0-161 ;  [4]  0068;  [6]  0023;  [6]  0006.  See  table  1  for  more  information. 
The  swell  ratio  a(oo)/a(0)  for  zero  surface  tension,  gravity  and  Reynolds  number  is  also  indicated 
(a( oo)/o,  =■=  1-113)  (see  Reddy  A  Tanner  1978;  Trogdon  A  Joseph  1980,  1981). 


this  problem  there  is  an  important  exit  region,  of  just  a  few  nozzle  diameters,  in  which 
the  rearrangement  of  the  velocity  profile  and  the  jet  diameter  to  their  final  constant 
values  is  nearly  complete.  This  exit  region  is  clearly  evident  in  the  photographs  in 
figure  2.  We  want  to  exclude  this  small  complicated  region  near  the  exit  from  the 
rest  of  the  jet,  which  has  a  simpler  dynamics. 

For  the  eh  acite  experiments,  the  swell  near  the  pipe  exit  is  in  good  agreement  with 
Tanner’s  (1970)  prediction  when  injection  is  into  air,  but  is  about  30%  greater  when 
injection  is  into  fluids  other  than  air.  The  same  effect  is  observed  with  other  visco¬ 
elastic  fluids  tested  and  appears  to  increase  with  the  elastic  shear  compliance  of  the 
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Figure  4.  Jet  shape  for elvacite/DEM  injected  into  air,  water  and  glycerol  solution:  [1]  equation 
(4.6);  [2)  ep  —  Sp  =  1067  (into  air);  [3]  <M) 67  (into  water);  [4]  0005  (into  glycerol  solution).  See 
table  3  for  more  information.  We  observe  a  very  large  die  swell  in  the  experiments  labelled  [3], 
[4].  The  die-swell  ratios  for  Newtonian  fluid  (Reddy  ft  Tanner  1978;  Trogdon  ft  Joseph  1980,  1981) 
and  for  viscoelastic  fluid  (Tanner  1970)  are  also  indicated  above. 


fluid.  At  one  extreme,  when  e  =  0,  the  flow  a  few  diameters  downstream  is 
independent  of  r  and  x.  At  the  other  extreme,  for  jets  of  liquid  into  air  with  e  «  1 
the  initial  swelling  does  not  appear  to  reach  its  die-swell  value  a,  because  the  effects 
of  the  acceleration  of  gravity  in  increasing  U  and  decreasing  a  are  manifest  even  quite 
close  to  the  exit.  On  the  other  hand,  jets  of  liquid  into  liquid  are  strongly  influenced 
by  the  change  of  momentum  of  ambient  liquid  dragged  along  by  the  jet,  and  the 
velocities  are  slower  and  the  jet  radius  larger  than  in  the  case  of  injection  into  air. 


3.  Momentum  of  the  jet  and  entrained  fluid 

The  equations  governing  the  flow  in  the  jet  are 

U  =  tx  U(x,r)+tT  W(x,r)  (velocity),  (3.1a) 

T  =  — fl  +  S  (stress),  (3.16) 

8r(rW)+dx{rU)  =  0,  (3.2a) 

p(WdrU+UdxU)  =  dxTxx+  -rdr(rTrx)+pg,  (3.26) 

p(Wdr  w+udx  W)  =  drTrr+dxTrx+  -( Trr-Tee ).  (3.2c) 

r 


The  same  equations  hold  outside  the  jet.  We  use  a  tilde  (V,  T,  p  etc.)  to  designate 
variables  in  the  fluid  outside  the  jet.  At  the  free  surface  r  =  a(x)  we  require 
respectively  that  the  normal  component  of  velocity  vanish  and  that  jumps 
[]  =  ()  —  Q  in  the  velocity,  the  shear  stress  and  the  difference  between  jump  in  the 
normal  stress  and  the  surface-tension  force  (<r  =  coefficient  of  surface  tension)  all 
vanish.  When  expressed  in  coordinate  form  these  conditions  may  be  written  as 

[0]  =  0,  (3.3a) 

W-a'U  =  0,  (3.36) 
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ITrr + a'*Txx — 2a’Txt\  *J' 

/2  *  _ 
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(3.3  c) 

1  _1 1  nn'  —  v’  (3.3(f) 

1  +a  1  aa 

where  J  =  a/ ( 1  -fa'2}),  and  —pi  is  the  stress  of  the  ambient  fluid  on  the  jet.  Equation 
(3.3c)  may  be  used  to  simplify  (3.3d): 


We  now  introduce  the  head 


lTrr-a'Txr}+~  =  0. 
aa 

®  =  p-pgx+Po> 

&  =*  p-pgx+p 0, 


(3.4) 

(3.5  a) 
(?  Si 


where  p0  is  the  pressure  in  the  ambient  fluid  at  the  pipe  exit  at  x  =  x0,  and  the  -a 
stress  S  into  our  equations.  Then,  in  the  jet  we  have 


p(  W8r  U+  U8X  U)  =  8X(  -  0  +  Sxx)  +  -rdr(rSrx),  (  , 

p(Wdr  W+  U8X  W)  =  dr(-4>+Srr)  +  dxSrx+  i(<Srr-^),  (3.65) 

whilst  on  r  —  a(x)  [U]  =  0,  (3.7a) 

W-a’U  =  0,  (3.7  b) 

a'lSrr-Sxx]  +  (1  -a'*)  lSrxJ  =  0,  (3.7c) 

— ep0*+[— ®+$rr— a'$r;rJ+  =  0.  {3.1  d) 


To  determine  the  hydrodynamics  of  the  jet  far  downstream  we  shall  need  to  relate 
the  extra  stress  and  velocity.  Analysis  shows  that,  though  the  velocity  becomes  large 
ultimately  like  0  <  a  <  J,  strains  decay  at  least  like  x~P,  fi  >  0.  It  therefore  seems 
appropriate  to  use  the  small-strain-rate  expressions  S  =  pAj-f  otjAg-f  a2  Af  of  the 
fluid  of  second  grade.  Here  p  is  the  viscosity, 

Ax  =  VU+VUT,  A,  =  (U .  V)  Aj+Aj  .  VU-f  (At .  VU)T, 
and  a,  and  at  are  constants.  In  component  form  we  have 
Srr  =  2p8r  W+a1[2W3rr  W+2U8rx  W+ 4(5r  W)2 

+  2 (8r  U)*  +  28rUdxW]+at  [4 {8r  W)*  +  {8r  U+ 8X  W)*l  (3.9a) 

S„.*&  +  tl(,WSrK  +  1Ul^+  4^)+4a,^,  (M.) 

8XX  =  2p8xU+a1[2U8xxU+2WdrxU+4(8xU)* 

+  2 (8X  W)*  +  28xW8rU]+ut  [4 (8X  U)*  +  {8rU+ 8X  W)*l  (3.9c) 

Srx  =  fi(8r  U+8XW)  +  ax  ( W(8 „  U+  8rx  W)  +  U(8rx  U+ 8XX  W) 

+  38rWdxW+3dTU8xU+drUdrW+dxWdxU] 

+a,{2dr  W8X  W+28r  Udx  U+28r  W8rU+28xW8xU].  (3.9  d) 


In  |5  we  will  avoid  making  assumptions  about  the  relation  of  stress  to  deformation 
in  the  hope  of  obtaining  good  jet-shape  equations  for  unsteady  problems  and  other 
problems  in  which  the  assumption  of  a  second-order  fluid  is  not  justified. 
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Equations  (3  2),  (3.6),  (3.7)  and  (3.9)  govern  the  motion  at  each  point  of  the  jet. 
We  shall  work  with  area-averaged  equations  expressing  the  conservation  of  linear 
and  radial  momentum  of  the  jet  derived  (Joseph  1980)  by  integrating  (3.6)  across 
the  cross-section  of  the  jet,  using  (3.2)  and  (3.7).  We  get 

[pU2  +  <S>-Sxx]rdr  =  aa'[-epgx+&]  +<rJ'  +  a[Srx-a'Sxx],  (3.10) 

dx  Jq 

a'~  f  [pUW-Srx]rdr  =  aa'[epgx-$]-<rJ'  +  aa'[Srr-a'Srx]  +  a’  f  [<t-«w]dr. 
dx J  0  Jo 

(3.11) 

Area-averaged  momentum  equations  also  hold  outside  the  jet.  Here  we  get 

[pfr  +  1b-Sxx\rdr  =  -a[Srx+a'@-Sxx)},  (3.12) 

—  [pWtj—  STX\rdr  =  o[$  —  <?rr  +  a'*?,.*]  4-  J*  [ff>  —  8eg\dr.  (3.13) 

We  may  eliminate  boundary  terms  by  adding  (3.12)  to  (3.10)  and  (3.13)  to  (3.11). 
Thus 

£{JW-U,*  +  J  [pV*  +  Q>-gxx\rdr^  =  -epqxaa'  +  <rJ\  (3.14) 

^[rWU-Srx]rdr+  ^{pm-STX]rdr} 

=  epgxa  —  +  [°  [<1>  —  £<#] dr  +  j*  [4>  —  §gg]dr.  (3.15) 

®  Jo  •'a 

The  terms  involving  ep  —  g(p  — p)  >  0  drive  the  flow.  The  jet  falls  because  it  is  heavy. 
But  (3.14)  suggests  that  the  dynamics  of  the  jet  can  be  strongly  influenced  by  the 
change  of  momentum  of  entrained  liquid  outside  the  jet. 


4.  Torricelli  formulas 

We  are  going  to  obtain  asymptotic  formulas  for  the  jet  far  downstream  from  the 
point  x„  of  extrusion.  For  these  large  values  of  x,  x—x0  ~  x  no  matter  what  may  be 
the  (fixed)  value  x0.  We  get  uniqueness  of  asymptotic  formulas  up  to  arbitrary 
translations  of  the  origin  of  x.  This  type  of  lack  of  uniqueness  comes  from  the  fact 
that  4>(x,  r)  is  determined  by  the  global  problem,  and  is  only  determined  to  within 
a  constant  by  asymptotic  analysis.  In  fact  the  equations  are  invariant  under  the 
transformation 

pD(x,  r)  j  pP(x-  x0,  r)  +  epgx„ j 

We  are  going  to  show  that,  if  the  ambient  fluid  is  inviscid  and  dynamically  passive, 
§=$  =  0,  there  is  a  unique  asymptotic  solution  that  perturbs  the  Torricelli  limit. 

The  Torricelli  limit  arises  for  large  values  of  x  whenever  the  ambient  fluid  is  inviscid 
and  dynamically  passive.  To  specify  this  limit  it  is  convenient  (but  definitely  not 
neccessary)  to  consider  the  1  inviscid  ’  flow  of  the  viscous  jet.  In  such  a  flow ,  the  fluid 
is  viscous  but  the  viscous  terms  in  the  equations  are  inactive,  as  in  potential  flow 
of  a  viscous  fluid.  The  flow  in  the  falling  jet  is  not  inviscid  at  any  finite  x,  but  it  tends 


I 


Jets  into  liquid  under  gravity  455 

to  an  inviscid  flow  as  x  -*  oo.  Suppose  that  the  fluid  is  really  inviscid.  Then,  the  extra 
stress  S  vanishes,  and,  if  we  also  assume  that  er  =  0,  (3.10)  reduces  to 


J*  (p(72  +  ®)rdrj  =  —  epgxaa' . 


The  Bernoulli  equation  for  this  flow  (along  streamlines)  is 


where 


^>[£/*+  W*]-epgx+<t>  =  C, 
ra(x) 

Q  —  J  U(x,r)rdr 


is  independent  of  x,  and  |IF|  =  o\U\,  U(x,r)  -»  u(x)  as  x  -»  oo  and  a(x)  -*  0.  Hence  in 
the  limit 

Q  =  ia2u.  (4.3) 

C  is  a  pure  constant  which  may  be  absorbed  into  <l>  and  hence  put  to  zero  in  (4.2) 
without  losing  generality.  So  (4.1)  may  be  written  as 


and,  using  (4.3), 


£  j*  (|pw2 — epgx)  r  dr  J  =  —  epgxaa' , 

-'■SI 


It  follows  now  from  (4.3)  and  (4.5)  that 

u(x)  —  (2  egxft  (4.6) 

is  given  by  Torricelli’s  formula,  modified  for  buoyancy.  We  get  the  same  solutions 
in  §5  when  x— *  00,  even  for  non-Newtonian  viscous  fluids. 


5.  Perturbation  of  the  Torricelli  limit 

The  asymptotic  expansion  for  large  values  of  x  of  the  functions  U(x,r),  W(x,r), 
<t>(x,  r),  a(x)  describing  the  flow  in  the  falling  jet  follows  from  the  observation  that 
the  continual  acceleration  of  the  jet  due  to  gravity,  coupled  to  the  condition  that  the 
mass  flux  is  constant,  forces  the  radius  of  the  jet  to  contract  as  x  increases.  In  the 
limit,  for  large  values  of  x,  the  variation  of  U(x,  r)  across  the  cross-section  of  the  jet 
must  tend  to  zero  so  that  the  leading  terms  are  independent  of  r.  In  fact,  it  is 
apparently  possible  to  compute  a  unique  asymptotic  solution  of  (3.2),  (3.7)  with 
5  =  $  =  0  and  (3.9)  in  the  form 

¥(a;,r)  =  Ir*»0Jn(x)=  I  r*»  X  *}>>  ar*1  -  CxM, 

n - 1  i-l 

o(x)  =  X  anx~ln,  C=feg)l,  o,  =  ,  (5.1) 

n  - 1  \eg/ 

where  1 

[U(x,r),  ( W^(x,r)]  =  ~[-dr'¥(x,r)tdx'¥(x,r)]. 

Equations  governing  the  gtn(x)  are  obtained  by  identifying  independent  powers  of 
r  in  the  equations  that  hold  in  the  interior  of  the  jet.  At  the  boundary  r  =  o(x)  we 
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use  the  expansions  for  a(x)  and  gtn(x).  Kaye  &  Vale  (1969)  also  used  the  representation 
(5.1)  in  their  study  of  Newtonian  jets  into  air,  but  their  actual  calculation  omits  the 
effect  of  viscosity.  We  have  carried  out  the  identifications  using  (5. 1 )  that  are  required 
to  establish  the  asymptotic  representations  (5. 19) — (5.23).  We  are  not  going  to  display 
the  identification,  which  is  tedious;  instead  we  use  the  integrated  equations  of 
momentum  (3.10)  and  (3.11)  to  determine  the  autonomous  equation  (5.18)  for  a(x). 
We  solve  this  equation  in  a  series  of  fractional  powers  of  x  and  then  trace  back  through 
the  derivation  to  obtain  the  other  flow  quantities. 

To  obtain  the  equation  governing  the  asymptotic  expansion  of  a(x)  for  large  x  we 
first  decompose  the  variables  appearing  in  (3.10)  and  (3.11)  as  follows: 


U(x,r) 

u(x) 

u(x) 

&(x,  r)~ 

-  W(x,  r)/r 

K(x) 

Jfl'(x) 

w(x,  r) 

<t>(x,r) 

y(x) 

y(x) 

f(x,r) 

Sgg(x,r) 

= 

T(x) 

6(x) 

+  r* 

0 

0 

+ 

t(x,r) 

6(x,r) 

8rr(x,  r) 

fl(x) 

0 

^(x,  r) 

Srx(x>r) 

r(x) 

0 

f(x,r) 

_  J(x)  _ 

_  a(x)  . 

0 

J(x,r)_ 

(5.2a) 
(5.26) 
(5.2c) 
(5.2  d) 
(5.2  e) 
(5.2/) 
(5.2  g) 
(5.2  h) 

The  functions  in  the  first  two  columns  are  of  leading  order  in  powers  of  x~i  and  the 
functions  in  the  last  column  are  of  higher  order;  P,  0, 12,  t  also  have  terms  proportional 
to  rl  but  the  ^-dependent  functions  of  proportionality  are  of  higher  order.  The 
assertions  can  be  verified  a  posteriori.  The  volume  flux 

Q  =  |°  U(x,r)rdr  >=]a*[u(x)+iatu(x)]+fl{ x)  (5.3) 

fo<*> 

/,(*)  =  Jo  fl(x,  r)rdr. 

Q *  =  ^a*[u*(x)+a,u(x)U(x)]+ft(x)  =  Jo*  U*{x,r)rdr+fs(x), 
where  (5.2)  and  (5.3)  define  /,(*)  and  /»(*)  uniquely.  We  compute 

lpU*+Q-SX3[]rdr  =  +  $a*[y(x)  +  ia*y(x)  ~  Hx)]  +/4(x), 

where  JA(x)  is  uniquely  defined  by  (5.4),  (5.2c)  and  (5 .2d).  Hence  (3.10)  becomes 

SxFoT  +  WtyW  + 1°*?  (*)  ~  r  (*)]}  =  —  epgxaa'  +  <ra'  +/»(x) ,  (5.6) 

where /s  =  J'(x)—  /4(x).  In  the  same  way,  we  compute  from  (5.2c, e) 

£  [<&(*,  r)-Sw]  dr  =  a[y(x)+$a*y(x)-0(a;)]  +/,(x),  (5.7) 


is  a  constant  and 
Moreover 


r< 


(5.4) 


(5.5) 


and  using  (5.2a,  6,/) 


J"  [pUW—Srx]rdr  «  —\pu(x) «'(x) a*(x)  +/,(*). 


(6.8) 
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Hence  (3.11)  becomes 

~i p4-(u'uaa)  =  epgxa  — cr  +  a[y+$a2y  — 0]+/8(x),  (5.9) 

ax 

where  /g(x)  =  -/7'(x) -J'/a. 

Turning  next  to  (3.7c)  with  <&  =  S  =  0,  we  find  that 

epgx  +  y  +  a*y  =  £l(x)  +  -  +/»(x),  (5. 10) 

a 

where/,  is  defined  uniquely  by  (3.7c)  and  (5.2).  From  (5.9)  and  (5.10)  we  can  solve 
for  y  and  find 

7  =  ^i/>^(a3««')+/io(^).  (511) 

where /10  are  all  the  higher-order  terms  that  arise  after  elimination.  Now,  using  (5.1 1), 
we  write  (5.9)  as 

d 

—  \pj^(a3uu’)  =  epgxa+a[y+fa*y - 0\— cr+fn.  (5.12) 

Using  (5.12),  we  eliminate  y  +  $a2y  in  (5.6),  and,  after  simplification,  find  that 

Ix{^aT  +a2(0~r)- =  ivga‘+im'+f !*•  (5I3) 

Using  (5.2)  and  (5.9),  we  find  that 

6— r  =  —  3/iu'— Za^uu’)’  —  3agit'*+/13.  (5.14) 

Hence 

+  «!<««')'  +  <**(«'*)] «*}'- &[<*(«»'«*)]'  =  \epga2  +  \(Ta'  +}u 

(5.15) 

20®  * 

Finally  we  write — - u(x)  —  — - — \a%u+f16,  (5.16) 

a* 

where  U  may  be  determined  from  the  shear-stress  condition  (3.3 b)  written  as 

2aw-|aw'-3aV+/,*  =  0.  (5.17) 

Since  the  leading  terms  in  U  and  a  are  of  order  xi  and  x~i  respectively,  H  is  of  order 
x~t.  We  can  therefore  use  it(x)  =  2 Q/a*  in  (5.15)  without  affecting  the  final  result.  It 
then  follows  from  (5.15)  that 


— 4 pQ*a'  -  \epgai  +  Q(iQ[a"a* —a'*a]  —  48a,  Q*  F — —  2  ^r-1 

La  a*  J 

+  12a1Q*Fow'-13—  + 20  + ipQ*(a'am  -  5aV  +  4a'3)  =  fcra»a'+/„(x).  (5.18) 

We  now  set  all  of  the  ft{x)  =  0,  i  =  1,  2,  17.  We  solve  the  resulting  equations 

as  asymptotic  power  series  in  x“i.  Then  we  can  verify  that  every  /<(x)  is  of  lower  order 
than  the  last  term  retained.  All  the  unknowns  may  be  generated  from  (5.18),  and  (5.1) 
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has  a  unique  solution  in  powers  of  x~i  with  leading  term  given  by  (4.5).  This  solution 
is  exhibited  as  (5.19) 

_R1Rt  3(^J  — f?s)  Rt 

°(X)  =  ^  +  T  +  ~2R^r 

2R\~^Rt Rt  R\  Ri  +  SRtRt/Rt  „  _ 

+  — M3 - +  32 3 - ? - +0<*1’  ,<U9“) 

_(2Q*\i  p  _  <r  p-AQ  r  _90|p  D  _  81  «1  » 

V  eg )  4epg  epg  ‘16 p  1  H2p 

The  coefficients  R1,...,Ri  are  positive  and  Rt  is  probably  positive.  The  function  u(x) 
and  the  leading  term  in  &(x)  can  now  be  computed  using  (5.17)  and  (5.16).  This  gives 

/  4  f.  2i?g  3RS  2  Rt  $RS—1R\  4J?2ri?6  ^  — §vl 

(5.196) 

u(x]  =  ^x-i+^x-i). 

Returning  now  to  (5.10)  and  (5.11)  we  compute 


(5.20) 


y  =  - 


w+4«*|j*>+4w(f!  -  ‘  §j)*4+  jL(w^  |  +  %>|)*- 

+  4epg[Rt(j*  W*H  +  0(*_i).  (5.21) 


7  =  r*[-&epgx-'  +  0(x  *)]. 


(5.22) 


The  error  terms  in  (5.2)  are 


'4‘ 

V*' 

vb 

x~* 

$ 

X~i 

f 

=  0 

x~^ 

6 

x~t 

ft 

x-\ 

A 

T 

arH 

.J. 

_x-l_ 

Following  Clarke  (1968),  we  note  that,  when  the  distance  y  =  x— x9  from  the  point 
y  =  0  of  extrusion  is  introduced  into  (5.19),  we  get 

a'F+7~ir+0(!ri) 

asymptotically,  for  large  y.  This  shows  that  the  asymptotic  solution  is  uniquely 
determined  only  up  to  the  first  two  terms  which  contain  the  principal  effects  of  inertia 
and  surface  tension.  Other  effects,  like  viscosity  and  extensional  stresses,  depend  on 
the  undetermined  constant  x0.  In  the  case  of  jets  of  liquid  into  air  there  is  good 
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agreement  between  experiments  and  the  y~i  law  for  values  of  y  some  pipe  diameters 
downstream  of  the  pipe  exit  at  y  =  0.  This  good  agreement  is  evident  in  the  graph 
labelled  [  1  ]  in  figures  3  and  4  and  in  some  results  published  by  Scriven  &  Pigford  ( 1959) 
and  Kay  &  Vale  (1969). 

The  assumption  $  =  §  =  0  of  this  section  is  not  good  when  e  <  1.  In  this  case  the 
rate  of  change  of  the  momentum  of  the  ambient  liquid  entrained  by  the  moving  jet 
can  be  even  larger  than  the  change  in  the  momentum  of  the  jet  itself.  When  p  is  small 
(e  *  1),  however,  as  is  the  case  with  air,  the  momentum  of  the  entrained  fluid  is  small 
and  does  not  strongly  influence  the  observed  dynamics.  This  is  the  reason  why  the 
analysis  given  in  this  section  matches  experiments  on  jets  of  liquid  into  air  but  not 
into  other  liquids. 

Even  if  we  confine  our  discussion  to  the  case  of  jets  of  liquids  into  air,  it  is  necessary 
to  maintain  a  distinction  between  the  jet-shape  equation  (5.18)  and  the  seven-term 
asymptotic  (for  large  x)  solution  (5.19)  of  (5.18).  The  jet-shape  equation  and  the 
improved  jet-shape  equations  derived  in  §6  evidently  apply  in  intermediate  regions 
where  the  flow  is  not  yet  asymptotic.  This  was  shown  by  Kaye  &  Vale  (1969),  who 
studied  a  special  case  of  (5.18)  for  air  (p  «  Sp)  based  on  the  assumption  that  r-variations 
of  quantities  on  the  left  of  (5.2)  are  zero  and  =  ag  =  0.  Then  they  integrated  the 
resulting  equation,  with  <r  =  0,  subject  to  experimentally  determined  initial  conditions 
for  silicone  oil  into  air.  They  got  very  good  agreements  for  the  jet  shape  in  regions 
of  not  so  large  x  where  the  asymptotic  formulas  fail. 

6.  Area-averaged  equations 

In  §5  we  used  area-averaged  equations  to  derive  the  asymptotic  form  of  solutions 
that  perturb  the  Torricelli  limit.  Now  we  shall  examine  more  general  approximations 
which  follow  from  area-averaging  of  the  equations  of  momentum  and  the  assumption 
of  weak  r-dependence.  We  again  confine  our  attention  to  the  axisymmetric  case,  but 
the  flow  may  be  unsteady.  For  unsteady  flow  we  must  add  pdU/dt  and  pdW/dt  to 
the  left-hand  side  of  (3.25,  c)  respectively.  Moreover,  the  condition  (3.35),  which  says 
that  the  jet  surface  is  a  streamline,  is  replaced  with  the  kinematic  condition 

W  =  ^+a'U,  (6.1) 

where  the  prime  denotes  differentiation  with  respect  to  x.  In  the  unsteady  case  we 
derive  the  following  equations  expressing  the  area-averaged  balance  of  mass,  axial 
momentum  and  radial  momentum  respectively: 

-o-  (6-2) 

^Ftf  rUdr+\_j  (P^+<^-lS**)rrfr]  *  -pegxaa’  +  <rJ'+a[3„+a'^-Sxx)]r^m, 

(6.3) 

rWdr+[  j](PWU~Srs)' r*r\  ’ 

=pegxa+  f  (<b-SM)dr-  +o[-$+8rr-a'8rx) .  (6.4) 

Now  we  make  our  basic  approximation,  which  may  be  described,  loosely,  as 
applicable  to *  thin  ’  jets  or  filaments.  More  precisely,  we  shall  work  in  to-be-determined 
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regions  of  the  jet  where  variations  of  the  axial  velocity  U(x,r,t)  and  ‘pressure’ 
<D(x,  r,  t )  are  approximated  by  the  first  two  terms  of  the  expansion  in  powers  of  r* 
(see  ( 5.2 )).  It  will  be  enough  to  say  that 

U(x,r,t)  ~  U{x,t)  +  r*d(x,t),  (6.5a) 

<D(x,  r,  t)  ~  Q(x,  t)  +  r*$(x,  t),  (6.5 b) 


and  „l.at  S(U)  is  determined  by  U  alone.  In  this  case  we  may  reduce  (6.2)-(6.4)  and 
(3.76,  c)  into  five  nonlinear  equations  in  the  five  unknown  functions  U(x,t),  fi{x,t), 
<D(x,0,  $(*,0, 

The  first  step  in  the  reduction  is  the  determination  of  W  corresponding  to  (6.5 a) 


using  (3.2a): 


W(x,r,t)  =  -  \rV'(x,  t)-\r*8'(x,t). 


(6.6) 


Since  U  is  now  determined  in 


terms  of  U(x,t)  and  tf(x,t),  S[U]  is  also  determined 


in  terms  of  these  two  functions.  For  example  the  stresses  in  a  fluid  of  second  grade 
may  be  determined  from  (3.0).  8XX,  Srr  and  STX/r  are  polynomials  in  r*  with 


coefficients  depending  on  x  and  t.  We  may  write  in  analogy  to  (5.2)  that 


Tm(x,0 

Ir-  6mM 
0  0 

JnM 


(6.7) 


Combining  the  expressions  (6.5)-(8.7)  with  (6.2H8-4),  we  get 

|««+|W+ £#]'«<>,  (6.8) 

PjtWU+&0]+ 

-  —  pegxaa'  -f  orj' + a[8xr + a'($ — 8xx)]r  _  #>  (6.9) 

-P^wu'+fa'd'} 

-[pfr*UU')+^(fru+2U'd)+kpa'dd'- I~ ^3tw]  ' 

_  N  atm+l  „ 

=  pegxa±<l>a+)fl>a'-'L-m--{6m--^J'+a[-&+5rr-a'5zr\r_a.  (6.10) 
The  interface  conditions  (3.7  c,  d)  reduce  respectively  to 

N  N 

•'is^-r„)+(i-oi«^«Tli  - — Sgg) + ( i — «'•) Sgf,  (6.ii) 

•  • 

pegz-<t>~ut&+  «*'%.)  +  « -t  +  Srr-a'Sr,, 

•  «* 

In  the  case  of  a  Newtonian  fluid  N  ■  1  and 

8^  -  2/t^-  - -MW+Vfi'l 


(6.12) 


■■w  JT!’  *t*m  '■ 
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Sxx  =  2n^  =  2ti[U'  +  r*fr], 

SM  =  2pW  =  -2M[\U'+yfrl 

S*r  =  *  [l£  +  =  AW-i^+2^)-Jr*^]- 

It  is  evident  that  $  may  be  eliminated  from  this  system,  leaving  four  equations  for 
four  unknown  functions  U(x,t),  Tl(x,t),  <D(x,  t)  and  a(x,t).  We  want  to  know  the 
conditions  under  which  these  four  nonlinear  equations  can  give  a  good  description 
of  the  dynamics  of  the  jet. 

Special  cases  of  the  equations  for  jets  of  liquid  in  air  (<Jp  =  p,  $  =  8  =  0)  which 
arise  when  t)  and  &(x,  t)  are  put  to  zero  have  been  considered  by  various  authors 
(Matovich  Sc  Pearson  1969;  Pearson  Sc  Matovich  1969).  In  this  r-independent  case 
we  get  « 


ja*+\a*U\  =  0, 


p^-haiU-[\ax(pVt  +  <t>- r0)]'  =  —  pegxaa'  +  aJ' , 


(6.13) 


(6.14) 


8  T  /  UU'  XT'  J' 

~pdi^>Ul^\P~2 - TVJ  ~  Pcgxu+<ba-ad0-<T~ ,  (6.16) 


o/(n0-r0)+(i~o'*)oT0  =  o, 


—  (egx+<t>)+&n— aa'r0+  — 7  =  0, 

aa 


(6.16) 

(6.17) 


where  T 0,  Q,,  80  and  t0  are  determined  by  U(x,  t).  Of  course,  it  is  not  generally  possible 
to  solve  these  five  equations  for  the  three  unknown  functions  U,  and  a.  The  common 
practice  is  to  ignore  (6.16)  and  (6.17).  (In  fact,  the  first  four  terms  of  the  asymptotic 
expansions  (5.10)  and  (5.20)  correspond  to  a  solution  which  may  be  derived  from 
(6.13) — (6.15)  in  the  steady  case  and  automatically  satisfies  (6.16)  and  (6.17)  to  the 
order  considered). 

In  the  steady  case  (6.13)-(6.15)  reduce  to 


Q  =  |a*l7(x)  =  a  constant  independent  of  x, 
—  ro)]'  =■  —pegxaa'  +  arj', 


(6.18o) 

(6.18b) 


— [kz*(lpUU' — 7,)]'  »  pegxa+<t>a—d9a— 


(6.18c) 


Eliminating  $,  we  find  that 


[ka'(pU'+69-r.)Y-[a{Wy>UU'--T')Yy  «  .  (6.19) 


For  second-order  fluids,  satisfying  (3.9),  we  get 

0,-r,  -  -WjtU'+(*l+*t)U''+<ixUir)t 

r,  -  Hstf'+aq  [V’W- (717"]— a,  17' IT). 


(6.20) 

(6.21) 


rtS*S*S&nr'< . 
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i 

a  (mm) 


Figure  5.  Jet  shape  for  silicone  oil  1  into  air: - ,  equation  (4.6); - ,  numerical  integration  of 

(6.22)  with  the  right-hand  side  set  to  zero  and  taking  two  observed  values  a(y  —  1 )  and  a(y  —  26)  as 
prescribed ;  + ,  experimental  data. 


The  average  ‘  extension*  1  ’  stress  for  the  steady  jet  in  regions  far  downstream  is  given 
by  (6.20).  Equations  (6.18)-(6.21 )  may  be  combined  into  one  single  nonlinear  jet-shape 
equation  for  a(x).  For  Newtonian  fluids  we  get 


“Wj  +  ~  $  -ipega'-frJ  -  -Jtra(^)  - 


(6.22) 


The  nature  of  the  simplification  leading  from  the  r-dependent  equations  (6.8M6.12) 
to  (6.22)  are  such  as  to  leave  the  region  of  validity  of  (6.22)  obscure.  Some  authors 
(Kaye  ft  Vale  1969)  have  treated  (6.22)  for  the  special  case  in  which  the  right-hand 
side  is  negligible.  The  asymptotic  analysis  of  $5  shows  that  the  first  four  terms  in 
the  asymptotic  expansion  are  asymptotic  solutions  of  the  equation  arising  from  (6.22) 
when  the  right-hand  side  is  set  to  zero. 

There  are  also  some  numerical  results  which  suggest  that  (6.22)  is  a  valid  equation 
under  some  circumstances.  In  particular,  some  good  results  have  been  obtained  from 
(6.22)  in  two  special  cases  of  jets  of  liquid  into  air.  Kaye  ft  Vale  (1969)  integrated 

-W~f<W?(£ -£)-»««■  (8-23) 


numerically  for  experimentally  given  initial  conditions  o(x0)  and  o'(x0)  measured  on 
jets  of  silicone  oil  into  air.  They  got  good  agreement  between  the  numerical 
calculations  and  the  observed  jet  shape. 

Kaye  ft  Vale  (1969)  also  did  experiments  with  a  Newtonian  jet  of  low  viscosity. 
But  they  could  not  get  the  jet  shape  from  (6.23).  (They  should  have  included  the 
surface-tension  term  } J .) 

We  have  also  integrated  (6.23)  for  jets  of  silicone  oil  1  into  air  (ji  *  10  P, 
ep  =  p  ws  9-43,  Q  ^  16g/s,  <r  =  212  dyn/cm).  We  use  a  two-point  method  for 
numerical  integration,  taking  two  observed  values  oy_,  and  av_u  as  prescribed  (see 
figure  5).  We  get  good  agreement  between  the  experimental  shape  and  the  computed 
one  even  in  regions  in  which  the  asymptotic  solution  fails.  But  we  did  not  succeed  in 
getting  numerical  computations  based  on  (6.23)  to  agree  with  the  jets  of  liquid  into 
liquid  because  the  momentum  of  the  entrained  liquid  is  important  when  the  density 
of  the  ambient  liquid  is  not  small. 
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7.  Boundary  layers  on  the  jet 

To  assess  the  effects  of  viscous  drag  on  the  momentum  of  the  jet  we  now  consider 
the  dynamical  effects  of  the  ambient  fluid  in  a  boundary-layer  approximation.  We 
shall  suppose  that  the  viscosity  of  the  jet  is  much  larger  than  the  viscosity  of  the 
ambient  fluid,  as  is  true  in  our  experiments.  In  this  case  the  shear  of  the  ambient 
fluid  may  influence  the  momentum  of  the  jet,  but  it  does  not  induce  strong  variation 
of  velocity  across  the  cross-section  of  the  jet.  We  shall  further  suppose  that  the 
ambient  fluid  is  Newtonian  and  that  the  boundary  layer  is  small  in  the  sense  that 
the  pressure  across  it  may  be  neglected. 

fc  =  0,  (7.1) 

and  x-derivatives  are  smaller  than  r  derivatives.  Then 


„dC  n dO  ,/d*D  1 dO 


dx 


dr 


\  dr* 


dr)' 


(7.2a) 


S-rC+  j-rff  =  0, 

dx  dr 


(7.26) 


where  P  —  fi/p  is  the  kinematic  viscosity  of  the  ambient  fluid.  We  require  that  0  and 
IP  be  continuous  at  the  jet  surface  r  =  a(x)  and  that  a'(x)  be  small. 

Using  these  approximations  we  may  write 

dO 

a(Sxr-a'8xx)  (7.3) 


and  (3.10)  becomes 


(pU*+Q- 8xx)rdr  =  —  epgaa'+pa~  (x,a)+<rJ‘.  (7.4) 


It  now  follows  directly  from  (7.2)  and  the  conditions 

C,  as  r  -»  oo 

,dV  ( x,a )  _ 

dr  dx)a 


that 


ap- 


rdr. 


(7.5) 


Equation  (7.5)  says  that  the  Bhear  force  at  the  jet  surface  is  equal  to  the  change  in 
the  momentum  of  the  entrained  liquid. 

Now  we  seek  to  determine  the  velocity  and  diameter  of  the  jet  far  downstream 
where  r-variations  of  *  </>  and  U  —  w(x)  in  the  jet  are  negligible.  Then  to  leading 
order  (3.1 1 )  may  be  written  as 

epgxaa'— <rJ'+aa'y  *  0,  (7.6) 

and  (7.4)  becomes 

l[ip«V+  ^  -Ja*S„]  -  toga'+jtaQ;  (7.7) 


G-jiw*. 

When  fi  «  0,  we  have  the  case  studied  in  f  4,  and  the  leading  asymptotic  balance 
associated  with  (7.7)  is  ^ 

£paV  -  epga*. 


(7.8) 
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The  right-hand  side  of  (7.8)  can  be  interpreted  as  the  weight  of  the  jet  per  unit  length, 
and  the  left-hand  side  gives  the  change  of  momentum  per  unit  length. 

We  will  show  that,  when  fi  4=  0,  the  leading  balance  is 

fl 

\epgd1  +fia-^  =  0.  (7.9) 

Equation  (7.9)  expresses  the  balance  between  the  weight  of  the  jet  per  unit  length 
and  the  shear  tractions  on  the  jet  wall. 

We  shall  study  the  boundary-layer  equations  by  a  Karman-Pohlhausen  method. 
First  we  evaluate  (7.2)  on  r  =  a(z).  We  have 


identically  in  x.  Hence 


u(x)  =  0(x,a(x)) 

„  v  80  ,80 

"(*)=&+“s7' 


-so  so  r .  ,atr 1 ,  ,.,ec 

because  w—a'u  =0  at  r  =  a.  Now  we  construct  an  approximate 

0( X,  r)  =  u(x)  (1  -i))af(tl,x), 
f(y,x)  =  1  +  Kt(x)ii*, 


satisfying 


A(x)  =  8(x)—a(x) 

0(x,  8)  =  0,  0'(x,  8)  *=  0,  0(x,a)  =  u(x) 

8*0  1  80  _  uu' 

dr*  +  a  8r  v 


(7.10a) 

(7.106) 

(7.10c,  d) 

(7.11  a,  b,  c) 


(7.11  d) 


on  r  *  a(x).  The  conditions  (7.11o-c)  are  satisfied  by  (7.10),  and  (7.11)  holds  if 

«'A*  A 
Kt(x)  =  —  -1  +  -. 

*'  2v  a 


Asymptotically  we  find  that  Kt  =  0(xi )  and 

Kt(x)  *  A/o, 

where  the  neglected  terms,  computed  a  posteriori,  are  0(1)  and  O(xH). 
Using  (7.12)  we  find 


0(x,  r)  *  u(x)  ( 1  —  7)*  [  1  +  £  7*J , 
fiad-£(x,a)-iiaH(z){- f)«-4*g. 


(7.12) 


(7.13) 


(7.14) 


The  leading  asymptotic  balance  (7.9)  of  momentum  in  the  jet  may  now  be  written 


(7.16) 
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6p  (i/cm1) 

Figure  6.  Variation  of  jet  radius  with  ip  - ,  theoretical  result, 

+  ,  experimental  data  for  elvacite/DEM  at  y  —  25  cm ;  O,  experimental  data  for  silicone  oil  1  at 
y  —  25  cm.  In  these  plots  the  origin  is  taken  at  the  orifice.  The  data  show  a  slope  of  — j,  in  agreement 
with  theory.  The  theoretical  line  can  be  made  to  pass  through  the  experimental  points  by  shifting 
the  origin  of  x. 


The  balance  of  momentum  of  entrained  fluid  is  given  by  (7.5)  and  may  be  evaluated 
using  (7.13).  We  first  compute 


pfrrdr*\  p^rdr  =  ptt*[A*(*  +  I^l  +  I^5)  +  oA(J  +  1fcift  +  ^if,)].  (7.16) 

*  J  a 


To  leading  order  (A*A|  =  0(xb)  with  neglected  terms  O(a^),  a  posteriori)  this  is 

1260a1  1260 W  a*  31S\v»/  (7I7) 

Now,  using  (7.14)  and  (7.17),  (7.5)  may  be  written  as 


taking  a(oo)  *  0,  we  get 


dzl$l5a"\epg)  J  ' 


(7.18) 
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Figure  7.  Comparison  of  the  jet  shape  of  liquid  injected  into  air  with  Torricelli  limit:  (o)  Torricelli 
limit,  «(y)  (2Qt/eg)~^  =  y~i;  (6)  experimental  data  for  silicone  oil  1;  (e)  experimental  data  for 
elvacite/DEM.  The  three  plots  tend  to  coincide  as  y-*  oo. 


a(y) -  .  ~  ~  y *• 

(W)1 

[2HS1  experimental  data  for  silicone  oil  I  injected  into  alcohol  solution;  [2]  ep  -  ip  «  0-161 ;  [3] 
04)68;  [4]  0023;  [5]  0005;  [6],  [7]  experimental  data  for  elvacite/DEM  injected  into  water  and 
glycerol  solution  ;  [8]  ep  —  ip  —  0067 ;  [7]  0006.  In  this  figure  the  origin  for  the  theoretical  curve, 
which  is  arbitrary,  is  taken  at  the  orifice.  We  count  theory  and  experiments  as  being  in  agreement 
when  it  is  possible  to  achieve  agreement  as  y  -»  oo  by  changing  the  origin  ofy  in  [1],  The  theory  agrees 
with  experiments  [2H4]  (silicone  oil  1)  and  [6]  (elvacite/DEM).  When  ip  is  very  small  (ex  >riments 
[5]  and  [7])  asymptotic  regions  are  far  from  the  orifice  (eventually  the  accelerating  i**t  n  >-t  thin 
out).  For  [6]  and  [7],  ip  ”  0005,  the  jet  diameter  has  not  yet  entered  a  region  u  asy.  ~>totic 
dynamics  suitable  for  comparison  with  [I  ]. 


1 
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Fiovrk  9.  Variation  of  jet  radiuti  with  /?: - .  theoretical. 


/I  171p*(?A 


a 


+  ,  experimental  data  for  elvacite/DEM;  O,  experimental  data  for  silicone  oil  I.  Jet  radii  are 
measured  at  y  =  25  cm  (the  origin  of  x  is  taken  at  the  orifice).  The  data  show  a  slope  of  approxi¬ 
mately  J.  The  theoretical  line  can  be  made  to  pass  through  the  experimental  points  by  shifting  the 
origin  of  x. 


and  when  the  viscosity  /t  of  the  ambient  fluid  is  zero 


(7.19) 


The  x“*  decay  is  much  slower  than  the  x“i  decay,  which  must  prevail  at  large  x 
when  fl  =  0.  The  asymptotic  powers  exhibited  in  (7. 18)  are,  however,  an  approximation 
at  best  which  is  based  on  the  notion  that  the  dynamics  of  the  entrained  fluid  may 
be  treated  by  boundary-layer  analysis,  and  that  this  analysis  itself  is  adequately 
approximated  by  the  Karman-Pohlhausen  procedure. 


8.  Comparison  with  experiments 

The  following  points  of  comparison  between  theory  and  experiment  merit 
emphasis. 

(i)  The  experiment  for  jets  of  liquid  into  air  are  in  good  agreement  with  the 
Torricelli  limit.  The  nature  of  these  agreements  has  already  been  discussed  at  the  end 
of  |5. 
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(ii)  The  theory  predicts  that  the  jet  radius  scales  with 


This  J-power  singularity  in  the  density  difference  holds  uniformly  in  the  viscosity  p. 
even  for  /?  =  0.  Therefore  the  J-power  law  is  more  robust  than  other  features  of  the 
dynamics  since  it  holds  when  the  dynamics  are  dominated  either  by  the  changes  of 
momentum  of  the  jet  (as  in  jets  into  a  vacuum,  or  air)  or  by  the  entrained  liquid 
outside  the  jet.  The  J-power  law  is  in  good  agreement  with  our  experiment 
summarized  in  figure  6. 

(iii)  For  large  values  of  x  the  jet  radius  is  proportional  to 

(20  ^  ^  i)- 

The  a  =  J  results  are  verified  for  liquids  into  air,  and  a  =  ^  appears  in  the 
observations  of  liquids  into  liquids  with  e  -+  0.  The  rate  of  change  of  momentum  of 
the  jet  (given  by  the  first  term  of  (7.7),  which  by  itself  leads  to  a  =  J)  is  still  important 
at  the  largest  values  of  x  in  the  experiments.  The  shear-stress  term,  which  is  equal 
to  the  change  of  momentum  of  the  entrained  liquid,  becomes  more  and  more 
important  as  x  is  increased  for  fixed  e  or  as  e  is  decreased  at  a  fixed  x,  with  a(e) 
decreasing  monotonically,  apparently  from  J  to  In  reservation,  we  note  that  the 
value  a  =  ^  is  not  yet  evident  for  the  smallest  e  *  0  005  because  the  diameter  of 
the  jet  was  still  increasing  at  the  largest  x  for  which  the  jet  was  stable  (see  figures  7 
and  8). 

(iv)  The  jet  radius  appears  to  scale  with  the  J  power  of  viscosity  as  predicted  by 
(718)  (see  fig.  9). 

Experimental  results  for  the  silicone  oil  2  are  in  the  same  measure  of  agreement 
with  theory  as  those  already  displayed  for  silicone  oil  1  and  will  not  be  displayed  here. 
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Supplement  to  Section  7 


< 


Recall  the  jet  shape  equation  (eqs.  6.22  and  7.16) 
for  the  problem  in  which  the  velocity  is  assumed  to  be 
independent  of  r 


2  a  »  a  " 

-4p(T  ^3  +  6yQ(~- 

aJ 


(i) 


2  2 
r.  „2  9  ,  K-  K~  ,  2K7  K- 

'j_E4-  [A  (I0  +  T70  +  T260}  +  aA(5  +  TTuT  +  ^Tn)  1 

d 


105  630' 


4yQ 

Aa 


(2) 


where 


K 


2 


1 


We  shall  now  show  that  these  equations  seem  good  for 
approximations  when  a  <<  1*  Then 


K 


2 


and  the  boundary  layer  equation  (2)  can  be  integrated 
analytically  to  give 


<£>5  +  W(J)4  +  45(J)  3  +  575  (J)  2  =  1575^  7 

a  lb  a  a  ay 


(3) 


Equations  (1)  and  (3)  can  be  simultaneously  solved  as 
a  two  point  boundary  value  problem  using  values  obtained 
from  experiments  at  x  =  1  and  x  =  i  as  prescribed.  The 
results  are  shown  in  Figures  7.1  and  7.2. 

Note  that  for  silicone  oil,  the  results  of  numerical 
calculations  agree  better  with  experiments  for  smaller  5p. 
This  is  due  to  the  fact  that  for  smaller  6p,  the  assumptions 
that  a  ~  constant  and  a'  <<  1  are  more  realistic. 


Figure  7.2  Jet  shape  of  elvacite/DEM: - ,  numerical  integration  of  equati 

(6.2)  and  (7.16)  taking  two  observed  values  a(y*=l)  and  a(y-f)  as  prescribed 
+  experimental  data. 

[1]  6p  ■  .067  (into  water);  [2]  Sp  =  1.067  (into  air). 
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